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ABSTRACT 

We describe photometric recalibration of data obtained by the asteroid survey 
LINEAR. Although LINEAR was designed for astrometric discovery of moving 
objects, the dataset described here contains over 5 billion photometric measure- 
ments for about 25 million objects, mostly stars. We use SDSS data from the 
overlapping ~10,000 deg 2 of sky to recalibrate LINEAR photometry, and achieve 
errors of 0.03 mag for sources not limited by photon statistics, with errors of 0.2 
mag at r ~ 18. With its 200 observations per object on average, LINEAR data 
provide time domain information for the brightest 4 magnitudes of SDSS survey. 
At the same time, LINEAR extends the deepest similar wide-area variability sur- 
vey, the Northern Sky Variability Survey, by 3 mag. We briefly discuss the prop- 
erties of about 7,000 visually confirmed periodic variables, dominated by roughly 
equal fractions of RR Lyrae stars and eclipsing binary stars, and analyze their 
distribution in optical and infra-red color-color diagrams. The LINEAR dataset 
is publicly available from the SkyDOT website (http://skydot.lanl.gov). 

Subject headings: binaries: eclipsing - RR Lyrae variable - stars: variables: 
general - astronomical databases: catalogs, surveys 

1. Introduction 

Variability is an important phenomenon in astrophysical studies of structure and evolu- 
tion, both stellar, galactic, and extragalactic. Some variable stars, such as RR Lyrae stars, 
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are an excellent tool for studying the Galaxy. Being nearly standard candles (thus mak- 
ing distance determination relatively straightforward) and be ing intrinsically b right, they 



are a particularly suitable tracer of Galactic structure (e.g., iSesar et al.l 120071. and refer- 
ences therein). Similarly, eclipsing binaries can be used as distance ind icators ( IGuinan et al. 
19981 ). are excellent probes of stellar physics (IPrsa fc Zwitterl 120051 ). an d offer a unique 



method for measuring stellar mas ses and other parameters as well (e.g., 
Torres. Andersen &: Gimenezll2010l ). 



Andersen 1991 



Despite the importance of variability, the variable optical sky remains largely unexplored 
and poorly quantified, especially at the faint end (V > 15). To what degree different variable 
populations contribute to the overall variability, how they are distributed in magnitude and 
color, and what are the characteristic timescales and dominant mechanisms of variability, 
are just some of the questions that still remain to be answered. To address these questions, 
several contemporary projects aimed at regular monitoring of the optical sky were started. 
The four most prominent wide- area (more than a few thousand deg 2 ) variability surveys in 
terms of depth and cadence are as follows: 



ROTSE-I (jAkerlof et al.ll2000l ) monitored the entire observable sky twice a night from 
V = 10 to a limit of V = 15.5. The Northern Sky Variability Survey (NSVS; 



Wozniak et al.l 120041 ) is based on ROTSE-I data and contains light curves for about 14 
million objects, with 100-500 measurements per object collected over one year. Typi- 
cal photometric precision of this dataset is about 0.02 mag for sources not limited by 
photon statistics. 



The All Sky Automated Survey (AS AS; iPojmahski 1120021 ) monitors the entire southern 
and part of the northern sky (5 < +28°) to a limit of V = 14, including about 15 million 
stars. Typical photometric precision is about 0.05 mag. The third release of the ASAS 
Catalog of Variable Stars contains close to 50,000 variable stars. About 80% of these 
are new discoveries (not previously cataloged). 



The Lowell Observatory Near Earth Objects Survey Phase I (LONEOS-I; iMiceli et al. 
20081 ) provides photometric data for 1430 deg 2 of northern sky that has been imaged 
at least 28 times between 1998 and 2000. The LONEOS-I camera used no bandpass 
filter and reached a depth of R ~ 18.5. Typical photometric precision of this dataset 
is about 0.02 mag for sources not limited by photon statistics. 



The Palomar Transient Factory (PTF fLaw et al.ll2009l and iRau et al.l 120091 ) is an on- 
going wide-area, two-band (SDSS-g' and Mould- R filters), deep (R ~ 20.6, g' ~ 21.3) 
survey aimed at systematic exploration of the optical transient sky. As of Fall 2011, 
PTF has observed about 7300 deg 2 of northern sky at least 30 times in the Mould- R 



-3 - 



band (~ 1800 deg 2 of sky with more than 100 epochs). Typical photometric precision 
of this dataset is better than 0.01 for sources not limited by photon statistics. 



A com prehensive review of past and ongoing variability surveys can be found in iBecker et al. 

fiooj. 



With hundreds of observations per object, the NSVS and ASAS datasets represent 
unique resources. Although a number of massive synoptic surveys , such as Pan-STA RRS 
(IKaiser et al.l 120021 ) and Large Synoptic Survey Telescope (LSST; llvezic et al.l 120081 ). will 
significantly improve our knowledge of the variable faint sky, they will not obtain hundreds 
of observations per object any time soon. 

A downside of NSVS, ASAS, LONEOS, and PTF datasets is the lack of precise multi- 
color photometry, whi ch can provide val uable information about sources in addit ion to light 



curve characteristics (jCovey et al.l 120071 ). The Sloan Digital Sky Survey (SDSS; lYork et al. 



20001 ) has obtained five-band photometry for hundreds of millions of stars, but its temporal 
coverage is very poor. The only region on the sky with mor e than ten observa tions per source 
is the ~300 deg 2 large equatorial region called Stripe 82 ( jSesar et al.l 120071 . and references 
therein). At the same time, NSVS and ASAS are not deep enough to make efficient use of 
SDSS photometry. This goal to combine a large-area variability survey providing hundreds of 
observations with exquisite SDSS photometry was recently made reachable by the LINEAR 
survey. 

MIT Lincoln Laboratory (MITLL) has operated the Lincoln Near-Earth Asteroid Re- 
search (LINEAR) program since 1998 to search for asteroids. The focus of the program is to 
discover and track near-Earth asteroids (NEAs) larger than 1 km diameter, but LINEAR has 
also discovered and tracked many main-belt asteroids and smaller NEAs. Since late 2002, 
the LINEAR program has maintained an archive of all of its imagery data. This archive 
now contains approximately 6 million images of the sky, most of which are 5 megapixel im- 
ages covering 2 deg 2 . The LINEAR image archive contains a unique combination of depth 
(V < 18), coverage area, and cadence that is potentially useful in exploring a number of 
time-varying astronomical phenomena. 

Since LINEAR project was focused on astrometric survey, standard stars required for 
precise photometric calibration were often unavailable. Furthermore, a significant fraction of 
data was taken in non-photometric conditions. The main purpose of this paper is to describe 
photometric recalibration of LINEAR data using stars observed by SDSS as photometric 
standards. Companion papers are discussed in § 



m 



In § [2] we describe LINEAR survey in more detail, and discuss recalibration procedures 
[3j The recalibrated dataset is characterized in § H] and a preliminary analysis of about 
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7,000 confirmed periodic variables is presented in § [5j Our results are summarized in § [61 

2. LINEAR Survey 

2.1. Hardware Characteristics 

The LINEAR program operates two telescopes at the Experimental Test Site located 
within the US Army White Sands Missile Range in central New Mexico, at an altitude of 
1,506 m above sea level, and a latitude of 33. 8° N. The program uses two, essentially identical, 
equatorial-mounted telescopes of the GEODSS (Ground-based Electro-Optical Deep Space 
Surveillance) type, usually denoted LI and L2. Each telescope is a folde d prime-focus desig n 



with a primary mirror diameter of 1.016 m diameter, and //# of 2.15 (IStokes et al.ll2000i ) 



The two telescopes of the LINEAR project are each equipped with a CCD camera de- 
veloped at MITLL using CCDs devel oped and manufactured by MITLL for the US Air Force 



GEODSS system flBurke et al.lll998f ). These CCID-16 sensors are 5 megapixel (2560x1960), 
back-illuminated, frame-transfer CCDs with fast readout, low dark current, and low read 
noise. The CCID-16 sensors have 8 readout regions (hereafter referred to as cells) and have 
high quantum efficiency over a broad part of the visible and near-infrared spectrum (with 
overall solar-weighted quantum efficiency of 65% and peak efficiency of 96% at 620 nm, see 
Figured]). The cameras have no spectral filters. The combination of the CCID-16 camera 
and the 1-m GEODSS telescope produces an instantaneous field-of-view of 1.60° by 1.23° 
(~ 2 deg 2 ), with 2". 25 pixels. In Section [31 the combination of CCID-16 sensors and LI and 
L2 telescopes will be referred to as the LI and L2 setups. 

For a few months in 2003 (approximately August to November) the camera on the L2 
telescope was replaced with a very similar camera containing a smaller format sensor called 
the CCID-10. This smaller sensor is identical to the CCID-16s in terms of manufacture 
and processing, readout electronics and settings, spectral quantum efficiency and pixel size. 
The only difference between the smaller format camera and the larger format camera is the 
number of pixels, 1024x1024 for the small camera, and the number of readout regions (4 
cells). During the four months in 2003 that L2 was operating with the CCID-10 camera, 
the field-of-view was about 0.64° by 0.64° (~ 0.4 deg 2 ). In Section [3J the combination of 
CCID-10 sensor and L2 telescope will be referred to as the small setup. 
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2.2. Observing Conditions and Cadence 

The telescope site has reasonably dark skies, approaching \iy — 21 mag arcsec -2 on the 
best nights. However, the LINEAR program operates on nights with suboptimal conditions 
including haze, bright moonlight (to within 3 days of the full moon), scattered clouds, 
and airglow. LINEAR images are frequently taken under skies as bright as \iy = 16 mag 
arcsec -2 . The point spread function includes unknown contributions from atmospheric seeing 
and various system contributions (e.g., defocus and dome airflow) and is typically around 5" 
FWHM, resulting in slightly undersampled images. 

The image integration time is allowed to vary over the seasons so that a roughly constant 
number of search fields fills the available dark time. During the earlier days of the LINEAR 
program (roughly 1998 to 2006), the integration time was 10-12 seconds on long winter 
nights and 3-5 seconds on short summer nights. In recent years (roughly 2007 to 2009), 
the integration time was increased to 15-18 seconds in the winter and 8-13 seconds in the 
summer. Under these observing conditions, the magnitude for which the magnitude error is 
~ 0.2 mag (the 5 — a detection limit) is about 18 mag on average. 

As shown in Figure [2j the cadence of LINEAR observations (time between revisits of 
the same patch of sky) is fairly uniform, with a peak at ~ 15 minutes and a gap at ~ 8 
hours. A more detailed description of LINEAR observing strategy is given in Appendix |A] 



2.3. Photometry and Astrometry using SExtractor 

The LINEAR project has been archiving all of its image data since December 2002 
through the present. The data used for this work encompass all imagery collected from 
December 2002 through March 2008. After selecting the LINEAR imagery data that overlaps 
the main SDSS survey region (i.e. galactic latitudes greater than 30°) and SDSS stripe 82 
region (20 h 32 m < a J20 oo.o < 04^00™, -1.26° < 5 J20 oo.o < +1.26°, ~ 280 deg 2 )), the data 
comprises about 1.8 million images. The bias and flat-field corrections were applied to these 
images following standard methodology. 

The extraction of sources from LIN EAR images and the cr eation of source catalogs 



were done using the SExtractor software (IBertin fc Arnoutslll996l ; see the SExtractor users 
manual] for more details on the source detection and measurement process). The prelimi- 
nary astrometry for extracted sources was obtained by converting centroid pixel coordinates 



1 Available at https : //www. astromatic .net /pubs vn/ software/ sextractor/trunk/doc/ sextractor.pdf 
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into equatorial J2000.0 right ascension (R.A.) and declination (Dec) by application of the 
gnomonic projection and a-priori telescope pointing information. This projection does not 
correct for telescope pointing errors on the order of several pixels, telescope optical distor- 
tion on the order of a few pixels near the corners of the image, or differential refraction 
on the order of a fe w pixels at typical airmasses. When compared to the USNO-B catalog 



( jMonet et al.l 120031 ) astrometry, the average uncertainty in the preliminary astrometry is 
about ~ 6" in R.A. and Dec coordinates. Therefore, the astrometry reported by SExtractor 
must be recalibrated before matching extracted LINEAR sources against other catalogs. The 
recalibration of preliminary astrometry is described in the next section. 

The preliminary fixed-aperture photometry for each source is computed by summing 
the pixels within a 5 pixel (~ 11") radius circular aperture and subtracting the background 
estimate for each pixel computed in the background estimate step. The uncertainty on the 
computed flux is computed assuming Poisson statistics in detected electrons with a gain 
conversion factor of 2.5e~ /ADU. The pixels with counts above 17000 were not used as 
the CCDs begin to experience non-linear response above this level (corresponding to 42000 
electrons), and they saturate at full- well capacity of 75,000 electrons. 



3. Recalibration of LINEAR astrometry and photometry 

In this Section, we describe the astrometric and photometric recalibration of LINEAR 
data and quantify the quality of recalibrated data. The recalibration and data quality 
assessment are done on a field-by-field basis, where a field is defined as part of a LINEAR 
image observed by a single CCD cell (there are 4 cells in the small setup and 8 cells in the 
LI and L2 setups). The only exception to this is the astrometric recalibration, which is done 
on a image-by- image basis. However, the quality of astrometric recalibration is still assessed 
on a field-by-field basis. 



3.1. Astrometric recalibration 



We recalibrate the preliminary astrometry using the Astrometry.net software ( jLang et al 



20101 ). The Astrometry.net uses the x and y (hereafter, pixel) positions of sources extracted 
by SExtractor and looks for "skymarks" (sets of four stars) that corresp ond to pre-computed 
"skymarks" found in the "clean" USNO-B catalog ( IBarron et al.l 120081 ). Once a skymark is 
found the code obtains an initial astrometric solution (a transformation from the pixel co- 
ordinates into the J2000.0 Equatorial system) in the form of a third-degree polynomial in 
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x and ?/ coordinates. Using this initial astrometric solution the code looks for additional 
reference USNO-B stars in the image, and when it finds them it further refines the initial 
astrometric solution. 

The quality of recalibrated LINEAR astrometry is estimated by comparing the LIN- 
EAR positions with reference USNO-B positions (Figure [3]). We find that the recalibrated 
LINEAR single-epoch coordinates have an rms scatter of ~ 0.6" and are offset < 0.3" with 
respect to the USNO-B coordinates. These systematic offsets are independent of pixel posi- 
tions or the position of the image on the sky, and are subtracted from recalibrated LINEAR 
coordinates. Using a subset of data, we have found that the offsets decrease if the Astrome- 
try.net is allowed to search for even more reference USNO-B stars when refining the initial 
astrometric solution. 

To find fields with bad astrometry, we calculate the median and rms scatter of astro- 
metric residuals (difference between LINEAR and USNO-B positions) in R.A. and Dec for 
each field, and tag fields as bad if the median or rms scatter in either coordinate (R.A. or 
Dec) is greater than 2". We find that about 7% of all fields have bad astrometry, with the 
majority of bad fields (~ 90%) originating from the L2 setup, and remove them from further 
processing. Visual inspection of images containing bad fields revealed that bad fields lack 
several pixel columns, which were most probably dropped due to cell readout problems. 

In order to recalibrate the LINEAR photometry, we positionally matched the LINEAR 



dataset to the Sloan Digital Sky Survey Data Release 7 (SDSS DR7; lAbazajian et al.ll2009f ) 
imaging catalog. The SDSS DR7 imaging catalog provides homogeneous and deep (r < 22.5) 
1-2% accurate photometry in five band-passes (u, g, r, i, and z) of more than 11,000 deg 2 
of the North Galactic Cap, and three stripes in the South Galactic Cap totaling 740 deg 2 . 
The SPS S positions are a ccurate to better than 0.1" per coordinate (rms) for sources with 



r < 20.5 (IPier et al.ll2003l ). and the morphological inform ation from the images allows reliable 



star-galaxy separation to r ~ 21.5 flLupton et al.ll2002l ). 



We match LINEAR sources to SDSS DR7 non-saturated, "primary" object^] using a 
3" radius, and find ~ 4.8 billion matched pairs (there are 24 million unique SDSS sources 
included in these pairs). There are ~ 1.6 billion LINEAR sources without a SDSS match 
("orphans"), and they include observations associated with stars saturated in SDSS (r < 13, 
~ 60% of orphans), image artifacts (e.g., due to charge bleeding and diffraction spikes), 
moving objects such as asteroids, and other transients. 



See http : //cas . sdss . org/dr6/en/help/docs/algorithm . asp?key=f lags 
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3.2. Photometric recalibration 

We model the recalibrated LINEAR magnitudes as 

munear = am instr + C (g - i) + Z(x, y), (1) 

where a measures non-linearity of SExtractor-supplied instrumental magnitudes mi nstr = 
—2.5log(counts), C is a color-dependent term that accounts for the per-exposure change in 
the effective central wavelength of the LINEAR bandpass, and Z(x,y) is our model for the 
field-specific magnitude zero-point dependent on the x — y position in the image, defined as 

Z(x, y) = Z Q + Z x x + Z 2 y + F(x, y), (2) 

where Z , Z\ y and Z 2 vary from field to field, and F(x,y) is the "super flat-field" (fixed for 
each observational setup). Note that the color- independent airmass term (gray extinction) 
is implicitly included in Z . 

To recalibrate the LINEAR photometry we need a catalog of photometric standard stars 
in the LINEAR photometric (unfiltered) system. Since such a catalog does not exist, we use 
the SDSS imaging catalog and utilize SDSS gri photometry to derive magnitudes in the 
LINEAR photometric system. 



3.2.1. LINEAR Magnitudes Synthesized from SDSS Photometry 

We define LINEAR magnitudes synthesized from the SDSS photometry as 

m sdss = r + f(g - i), (3) 

where f(g — i) is some function of g — i, and SDSS gri magnitudes are not corrected for the 
ISM extinction. 

We find f(g — i) iteratively by first adopting a = 1.0 in Equation [T], F(x,y) = in 
Equation [2, and f(g — i) = in Equation [3j For each field, we determine Z , Z\ and Z 2 by 
minimizing 

X S[(t?7 [inear ^sdss) I ' 0~ linear] 5 (4) 

where o~ii near is the magnitude error supplied by the SExtractor (errors in m s d ss are assumed 
to be much smaller than errors in mu near ). 

The LINEAR observations of calibration stars must be matched within 2" to isolated^!, 
unresolved, SDSS objects brighter than r = 17, with < g — i < 2.6. Having Z(x, y) for 



3 SDSS CHILD flag set to 
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each field, we calculate mi nstr + Z(x,y) — r residuals for all calibration stars, bin the values 
in g — i bins, and plot the median residuals in Figure S] for the small, LI, and L2 setup. We 
plot medians for each setup separately because each setup represents a separate optical, and 
therefore photometric, system. 

Even though the three setups are different, their residuals depend similarly on g — i 
(rms scatter around the mean is 0.01 mag), as indicated by medians plotted in Figure HI We 
therefore fit a fourth-degree polynomial to all medians, and obtain 

f( g -i) = 0.0574 + 0.004(# - i) - 0.056(# - if + 0.052(# - if - 0.0262(# - if. (5) 

The transformation from SDSS gri to LINEAR magnitudes, as defined by Equation [3j is 
then 

m sdss = r + 0.0574 + 0.004(# - i) - 0.056(# - if + 0.052(# - if - 0.0262(# - if. (6) 

This equation effectively turns the SDSS imaging catalog into a catalog of LINEAR photo- 
metric standard stars. 

3.2.2. Super Flat-fields and Non-linearity of Instrumental Magnitudes 

We find super flat-field, F(x,y), iteratively by first adopting a = 1.0 in Equation [T] and 
F(x,y) = in Equation [2j For each field, we determine Z , Zi, Z 2 , and C by minimizing 
Equation H] using calibration stars brighter than m S( i ss = 17. The super flat-field correction, 
F(x, y), is then obtained by finding the median of mu near — m s d ss residuals binned in x — y 
pixels, and is shown as a color-coded map in Figure [5] (left) for each observational setup. 
Figure E] (left) shows the residuals after the subtraction of Figure [5] (left) maps. On average, 
the median residuals are now smaller than 0.005 mag, compared to ~ 0.02 — 0.03 mag before 
the correction was made. 

To verify our assumption of a = 1.0 (the linearity of instrumental magnitudes), we 
adopt a = 1.0 in Equation [TJ and determine Z , Z\, Z 2l and C by minimizing Equation @] 
using calibration stars brighter than m s d ss = 17. We bin mu near — m s d ss residuals in m s d ss 
bins, and plot the medians as a function of m s d ss m Figure [71 A strong dependence of 
residuals on m s d ss would indicate non-linear (a ^ 1.0) behavior of instrumental magnitudes. 
With the exception of a small non-linearity (< 0.02 mag) in the LI setup for m s d ss < 15, the 
residuals do not depend on magnitude indicating that the a = 1.0 assumption is valid. 

To estimate the magnitude of correction introduced by including the color-dependent 
term C (accounts for the change in LINEAR bandpass due to varying airmass) in Equation [TJ 
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we measure the rms scatter of C{g — i) values obtained from all calibration stars and fields. 
We find the rms scatter to be ~ 0.02 mag indicating that by correcting for this term the 
systematic uncertainty in LINEAR photometry is reduced by ~ 0.02 mag. On the other 
hand, the inclusion of the Z term is much more important as it reduces the systematic 
uncertainty in LINEAR photometry by ~ 0.3 mag (the rms scatter of Zq values is ~ 0.3 
mag). 

Once the color-dependent term, C, and the zero-point, Z(x,y), are obtained for a field, 
the LINEAR photometry for that field is recalibrated using Equation [TJ For "orphan" 
sources (LINEAR sources not associated with SDSS objects) the color-dependent term is 
set to C = 0, as g — % colors are not available for these sources. Therefore, the systematic 
uncertainty in recalibrated photometry is slightly higher for orphans (by ~ 0.02 mag). 



3.2.3. Checks on Recalibrated Photometry 

To quantify the quality of recalibration on a field-to-field basis, we calculate mu near — 
m s dss residuals for calibration stars, and find their median and rms scatter (determined from 
the interquartile range) for each field. The median of (mu near — m s d S s) residuals per field 
is then an estimate of the zero-point error. The distribution of these median offsets, shown 
in Figure El is about 0.01 mag wide and demonstrates the overall quality of recalibrated 
photometry. About 10% of fields are of low-quality (rms scatter in m« near — m s d ss > 0.1 
mag), usually due to highly variable cloud coverage. Sources detected in such fields are 
removed from further consideration. 

Another useful statistic is the rms scatter of \ values, {rn Unear — m sdss ) / au neari in m sdss 
bins. This rms scatter should be equal to 1 and should not depend on m s d ss if the SExtractor- 
supplied magnitude errors (au near ) are correct Gaussian errors. As shown in Figure [9j the 
rms scatter is greater than 1 indicating underestimated o~u near errors by up to a factor of 
~ 2. 

We correct underestimated magnitude errors by multiplying them by 1.4, 1.85, and 1.65 
for the small, LI, and L2 setup, respectively. In case of the small setup, the factor of 1.4 
is only the average value as SExtractor-supplied magnitude errors for that setup seem to 
depend on magnitude. Therefore, the magnitude errors for the small setup may be slightly 
overestimated or underestimated depending on the brightness of the source. This is not a 
major issue as only a small fraction of all LINEAR sources (~ 0.1%) originate from the small 
setup. 

The reason for this underestimation of SExtractor-supplied magnitude errors is not 
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clear. Non-inclusion of systematic errors (e.g., due to flat-fielding) in au near cannot explain 
this behavior for LI and L2 setups. If additive systematic errors were the reason for this 
underestimation, then (mn near —m s dss) / ^linear should decrease with magnitude because au near 
is expected to increase with magnitude (see Figure [T51) . Yet, no such decrease is seen in 
Figure [9] for LI and L2 setups. 
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4. Characterization of the Recalibrated LINEAR Dataset 

There are about 5 billion sources in the recalibrated LINEAR dataset. These sources 



were grouped into clusters using an implementation of the OPTICS algorithm (lAnkerst et al. 



19991 ) and a clustering radius of 2". There are about 25 million clusters with more than 15 
observations. SDSS provides morphological classification for 24 million clusters (7 million 
are resolved and 17 million are unresolved by SDSS), and for 1 million of these clusters there 
is no SDSS classification (they are slightly outside the SDSS footprint). Clusters with 15 or 
more sources are hereafter labeled as LINEAR objects or time-series. The median number of 
epochs (observations) per object as a function position and magnitude is shown in Figures [TD] 
and[TTJ The average number of epochs for objects brighter than 18 mag within ±10° of the 
ecliptic plane is ~ 460, and ~ 200 elsewhere. 

Clusters with less than 15 sources are not tagged as objects. We have statistically 
analyzed properties of 385 million sources (~ 8% of all LINEAR sources) assigned to such 
clusters and have found that most of them are very faint (mu near > 19.5), have very uncertain 
magnitudes {au near > 0.4 mag), and have high ellipticity (suggestive of cosmic rays and image 
artifacts). Such clusters and their observations are not included in our database. This also 
means that the LINEAR database does not include observations of fast moving objects (i.e. 
objects moving faster than ~ 1" day -1 ). 

In order to enable discovery of variable objects and to further characterize the LINEAR 
dataset, we have computed various low-order statistics for time-series, such as the median 
recalibrated magnitude ((mu near )), rms scatter (a), \ 2 P er degree of freedom (Xpdf)i skewness 
(7!), and kurtosis (72). The cumulative distributions of Xpdf f° r LINEAR objects with more 
than 30 observations are shown in Figure [T2J 

The median of the x\df distribution for bright LINEAR objects (14.5 < (mi inear ) < 17) 
is about 1, indicating that, on average, the uncertainties in LINEAR magnitudes are well 
determined. For fainter LINEAR objects (17 < (m linear ) < 18), the median of the Xpdf 
distribution is slightly lower (~ 0.7) indicating slightly overestimated photometric errors. 

The high-end tail (Xpdf > 3) of the cumulative x^df distribution suggests that ~ 5% 
of bright LINEAR objects exhibit deviations from the mean that are inconsistent with the 
Gaussian distribution of errors. These objects could be variable or could have unreliable 
(non-Gaussian) errors. While it is not impossible that some fractio n of LINEAR obj ects 



are variable (~ 1% of sources brighter than r = 20 are variable, see ISesar et al.l 120061 and 
references therein), a fraction of ~ 5% seems unlikely. A more likely explanation is spurious 
variability induced by non-Gaussian errors, caused by blending of sources and unaccounted 
instrumental defects (e.g., bad pixels). 
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T o estimate the dep endence of the median photometric error (£) on magnitude, we 



follow ISesar et al.l (120071 ) and calculate median a values in (mn near ) bins. The dependence 



of median a on (m[ inear ) , shown in Figure [£31 can be modeled as a fourth-degree polynomial 

£0) = 122.995 - 30.346a; + 2.8183a; 2 - 0.116893a; 3 + 0.0018295a; 4 , (7) 

where x = mi inear . The smallest photometric errors are at m linear ~ 15 (~ 0.03 mag) and 
are limited by systematic errors in calibration. The increase in photometric error at brighter 
magnitudes is probably due to saturation issues in LINEAR. The photometric error starts 
to increase rapidly towards fainter magnitudes due to photon noise, and reaches ~ 0.2 mag 
at mi inear ~ 18 (the adopted faint limit). 

In addition to statistics calculated on light curves, we also calculated time-averaged 
positions (median and mean right ascension and declination). To estimate the reliability of 
time-averaged LINEAR positions, we compared them to SDSS positions in ~ 1 deg 2 pixels 
on the sky and calculated the median difference in positions for each pixel. We find that, 
on average, the time-averaged LINEAR positions are within 0.3" of SDSS positions, with 
no systematic trends across the sky. The largest median difference (for a ~ 1 deg 2 pixel) 
is ~ 1".5. While this precision is sufficient for most applications, we still recommend that 
SDSS positions are used whenever possible to avoid patches where time-averaged LINEAR 
positions may be less reliable. 

Since SDSS is about 4 mag deeper than LINEAR, it can be used to estimate the com- 
pleteness of the LINEAR object catalog (i.e., LINEAR object catalog is a subset of the SDSS 
object catalog). To estimate the completeness of the LINEAR object catalog, in Figure [T~4l 
we plot the fraction of SDSS objects matched to LINEAR objects as a function of m S( i ss 
(LINEAR magnitude synthesized from SDSS photometry, see Equation [6]). For this purpose 
we use SDSS and LINEAR objects with < g — i < 2.5 (color range where Equation [6] 
is valid), 170° < RA. < 220°, and 15° < Dec < 65°. As shown in Figure EJ the LIN- 
EAR object catalog is more than 90% complete for m s d ss < 19. This estimate may be a bit 
misleading as the criterion for object creation is having at least 15 detections. Even faint 
objects (18 < m sdss < 19.5) are likely to be detected at least 15 times out of an average of 
200 observations (see Figure ITT]) . if the observing conditions are better than average (e.g, 
due to better seeing and darker skies). 
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5. Preliminary Analysis of Variable Objects Selected from the LINEAR 

Dataset 

We have examined a small subset of objects with Xpdf > 3 (~5,000) and identified a 
subset of periodically variable stars. We confirmed expectation from the preceding Section 
that only a few percent of the x%df > 3 subsample are convincing cases of periodic variability. 
As an illustration of the quality of LINEAR light curves at the bright {(rnu near ) ~ 15) and 
faint end {{mu near ) ~ 17), we show their period-folded light curves in Figure [151 The main 
strength of the LINEAR survey is easily discernible in these plots: even at the faint end, the 
light curve features are well defined thanks to dense LINEAR sampling. 

To exploit the full potential of the LINEAR photometric database, an automated se- 
lection and classification of variable objects is desirable. This desideratum is not easy to 
accomplish because the fraction of spurious objects with Xpdf > 3 is an order of magnitude 
larger than the fraction of truly variable objects. In order to enable the deployment of au- 
tomated selection and classification methods, we have undertaken an extensive program of 
visual light curve classification. Details of this program are described in a companion paper 
(Palaversa, L. et al. 2011, in preparation, hereafter Paper II) and here we only provide a brief 
description of the variable candidate selection and their distribution in various color-color 
diagrams. 

In the first step, we define a flux-limited sample of objects with 14.5 < (mu near ) < 17 and 
select about 200,000 light curves with x\df > 3 and a > 0.1 mag (cr is the rms scatter in the 
light curve). We limit classification to objects exhibiting periodic variability and use phased 
light curve s for visual in spection. The periods are determined using the Supersmoother 



algorithm flReimannlll994h . In total, about 7,000 LINEAR objects (3.5% out of 200,000) were 
visually confirmed as periodic variables, with about two thirds classified as RR Lyrae type 
ah stars and eclipsing binary stars of all types. The remaining 1/3 include SX Phoenicis and 
Delta Scuti stars, long-period variables, as well as possible RR Lyrae type ab and eclipsing 
binary stars with less certain classification. The distributions of classified objects in the 
SDSS u — g vs. g — r color-color diagram is shown in Figure [TBJ 

The differences between the distributions of all sources and those of the confirmed 
variable subsample demonstrate that the latter are robustly classified (i.e., the classification 
is not random). The most obvious difference is a much higher fraction of RR Lyrae stars 
in the classified sample (u — g ~ 1.15, g — r < 0.3). The eclipsing binary systems of all 
types dominate the visually classified sample in the main stellar locus region. There seems 
to be a hint for the median period of eclipsing binaries becoming longer for bluer stars. 
The distribution of eclipsing stars is offset from the distribution of all stars in the direction 
perpendicular to the stellar locus (offset to upper left). There is also a lack of variable 
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(eclipsing) systems of the late K and early M spectral type. These results are discussed in 
more detail in Paper II. 



6. Summary 

In this work we have described astrometric and photometric recalibration of data ob- 
tained by the asteroid survey LINEAR. The public access to the recalibrated LINEAR 



dataset will be provided through the SkyDOT Web site (http://skydot.lanl.gov). A 
description of various cataloged parameters that are available in this database is given in 
Appendix [B] 

The preliminary astrometry produced by the survey has been recalibrated using the 
Astrometry.net software and USNO-B catalog as the reference astrometric catalog. The 
precision (rms scatter) of single-epoch recalibrated LINEAR coordinates is ~ 0".6 with 
respect to the reference (USNO-B) astrometric catalog. The precision of time-averaged 
LINEAR coordinates is ~ 0".3 with respect to SDSS positions. 

In order to recalibrate the LINEAR photometry, we positionally matched the LIN- 
EAR dataset to the SDSS DR7 imaging catalog and used SDSS as a catalog of photomet- 
ric standard stars. LINEAR instrumental magnitudes were corrected for changes in the 
LINEAR bandpass (due to varying airmass) and for zero-point variations across the CCD 
(due to imperfect flat-fielding). These corrected LINEAR magnitudes were then offset on a 
per-LINEAR field basis to match the SDSS magnitude zero-point system. In general, the 
recalibrated LINEAR magnitudes (mi inear ) behave linearly, with the exception of a small 
non-linearity (< 0.02 mag) in the LI setup at magnitudes brighter than 15 mag. 

We find the average uncertainty in recalibrated LINEAR magnitudes to be ~ 0.03 mag 
at the bright end (mn near ~ 15, the systematic uncertainty) and ~ 0.2 mag at the faint end 
{miinear ~ 18, the adopted faint limit for the survey). The SExtractor-supplied magnitude 
errors were found to be underestimated by 40% — 85%, independent of magnitude, and were 
subsequently corrected for this behavior using multiplicative factors. 

There are about 5 billion photometric measurements in the recalibrated LINEAR dataset. 
The sources detected in each observation were grouped into clusters (time-series or objects) 
using a clustering radius of 2". There are about 25 million clusters with at least 15 observa- 
tions (7 million resolved and 17 million unresolved by SDSS, the rest without morphological 
classification). The median number of epochs (observations) per object brighter than 18 
mag and within ±10° of the ecliptic plane is ~ 460, and ~ 200 elsewhere. The completeness 
of the LINEAR object catalog (relative to the SDSS imaging catalog) was estimated to be 
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better than 90% for sources brighter than 19 mag. 

LINEAR dataset represents a major new resource for studying the variability of faint 
optical sources. With an average of 200 observations per object, LINEAR data provide time 
domain information for the brightest 4 magnitudes of the SDSS survey. At the same time, 
LINEAR significantly extends the deepest similar wide-area variability survey, the Northern 
Sky Variability Survey, by 3 mag. 

In order to help automated discovery of variable objects, we have computed various low- 
order statistics for time-series (median recalibrated magnitude, rms scatter, \ 2 P er degree of 
freedom, skewness, and kurtosis). The Xpdf distribution for bright LINEAR objects (brighter 
than 17 mag) is centered at ~ 1 (indicating that, on average, the uncertainties in LINEAR 
magnitudes are well determined), and has a long high-end (Xpdf > 3) tail. About 5% of 
LINEAR objects have Xpdf > 3; in a sample with Gaussian errors, less than 0.01% of objects 
are expected to have x\df > 3. This high-end tail is populated by truly variable objects, 
but is most likely dominated by spuriously variable objects with unreliable (non-Gaussian) 
errors (caused by, for example, unaccounted instrumental defects such as bad pixels or by 
blended sources). 

Since the fraction of spurious variable objects is at least a few times greater than the 
fraction of truly variable objects, we have undertaken an extensive program of visual light 
curve classification in order to enable future deployment of automated methods. Details of 
this program are described in a companion paper (Palaversa, L. et al. 2010, in preparation). 
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and recommendations are those of the authors and are not necessarily endorsed by the United 
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A. LINEAR Observing Strategy 

To search for asteroids, the LINEAR program attempts to cover the entire portion of 
sky available from the site each lunar month (lunation), with repeat searches within 10 to 15 
degrees of the ecliptic plane. The program schedules searches on both telescopes an average 
of 25 nights per lunation, taking off 4-5 nights around the full moon. After accounting for 
weather and equipment problems, about 20 observing nights are achieved per lunation with 
each telescope. To minimize stray light, the survey fields are kept about 20° away from the 
moon during the week before and after the full moon. 

To efficiently cover the available sky at minimal airmass, each telescope is assigned a 
strip of sky that is aligned with either the ecliptic or the equator and is typical 3-4 fields 
of view (4 — 5°) wide in declination or ecliptic latitude, and 150 fields of view (240°) long 
in right ascension or ecliptic longitude. The search box is divided into 5-7 regions each 
containing about 100 fields. The regions are searched starting with the westernmost region 
at the beginning of astronomical twilight in the evening, and ending with the easternmost 
region at the end of astronomical twilight in the morning. All of the fields in a single region 
are imaged in sequence, and then the sequence is repeated in that region 4 more times. The 
5 images of a given field are combined into a "frameset" which forms the input to a moving 
target detection algorithm that distinguishes asteroids from stars and noise based on their 
motion. This process of cycling through the fields in a region 5 times produces a temporal 
separation of 15 to 20 minutes between successive images of the same patch of sky. 

On subsequent nights, the entire search area for each telescope is generally stepped 
north or south by the height of the region in order to cover the entire sky, though gaps 
are occasionally unavoidable due to missed nights. Over the course of one lunation the 
entire available sky is generally searched once, and areas within about 10° of the ecliptic 
are searched two or three times. Because the monthly search regions are quite wide in right 
ascension, there is significant overlap in coverage area on adjacent lunar months. Therefore, a 
given patch of sky can be observed up to 8 months in a row before it passes into the daytime 
portion of the sky for 4 months. After a full year the search covers the entire sky north 
of declination —30° except for a small circle north of declination 82° where the telescope 
has poor settling ability. LINEAR's pattern of sky coverage produces repeat photometry 
data for stars on several timescales ranging from the 15-20 minute interval between images 
within a frameset, to a few days between repeat visits during one lunation, to the month- long 
timescale between lunar months, to the yearly. 
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B. LINEAR Database Schema 

The LINEAR database hosted at SkyDOT contains two tables, "Object" and "Source" . 
The "Source" table contains photometry and astrometry of individual LINEAR detections 
(sources), while the "Object" table references time-series obtained by clustering sources 
within 2". In addition, the "Object" table contains various low- level statistics calculated on 



time- series and is cross-matched to external catalogs such as SDSS and 2MASS (ISkrutskie et al 



20061 ). A detailed description of various parameters provided by these tables is given in Ta- 



bles [1] and [21 and a screenshot showing the SkyDOT interface is shown in Figure [T71 
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Table 1. LINEAR "Source" Table Schema 



Field 


Description 


sourcelD 


A unique number identifying this source 


ra 


Equatorial J2000.0 right ascension 


decl 


Equatorial J2000.0 declination 


MJD 


Modified Julian Date 


mag 


Recalibrated LINEAR magnitude 


magErr 


Uncertainty in the recalibrated LINEAR magnitude 


calib 


1 if the source was used for calibration, otherwise 


cllip 


"Ellipticity" parameter measured by SExtractor 


ximg 


Position of the source along CCD's x axis 


yimg 


Position of the source along CCD's y axis 


FWHM 


Full- Width- At-Half-Maximum measured by SExtractor 


flags" 


SExtractor internal flags 


objectID 


A number associating this source with an object in the "Object" table 



a See 



p. 28 in the SExtractor manual 



: //www . astromatic .net/pubsvn/sof tware/sextractor/trunk/doc/sextractor . pdf 
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Table 2. LINEAR "Object" Table Schema 



Field 


Description 




objcctID 


A unique number identifying this object 






Mean equatorial J 2000.0 right ascension 




IT 


Mean equatorial ,72000.0 declination 




ra_std 


Rms scatter in right ascension (in units of arcscc) 




dccl_std 


Rms scatter in declination (in units of arcscc) 




r a_m c d i an 


Median equatorial ,12000.0 right ascension 




d c cl_m cd i an 


Median equatorial J 2000. declination 




nptsTotal 


Number of detections 




nptsGood 


Number of detections with flags— — 




XXX ag _m can 


Mean recalibrated LINEAR, magnitude 




m. ag _w m. can 


Weighted recalibrated LINEAR magnitude {by magnitude errors) 




mag_iiicdian 


M^cdian recalibrated LINEAR magnitude 




mag_std 


Rms scatter of recalibrated LINEAR magnitudes 




mag_wstd 


Weighted rms scatter of recalibrated LINEAR magnitudes (by magnitude errors) 




mag_skc w 


Skcwncss of recalibrated LINEAR, magnitudes 




mag_kurtosis 


Kurtosis of recalibrated LINEAR, magnitudes 




mag_chi2dof 


per degree of freedom around the mean magnitude 




m ag_r ch i 2 d o f a 


R.obust x 2 P cr degree of freedom around the mean magnitude 




'< i — 5 

thcta_sdss 


Distance in arcscc from the SDSS source 




sdssID 


SDSS Data R.clcasc 7 objID parameter 




sdss_ra 


SDSS J2000.0 right ascension 




sdss_dcc 


SDSS J2000.0 declination 




objtype 


6 if unresolved in SDSS imaging (i.e., a point source), 3 if resolved 




rExt 


Extinction in the SDSS r-band 




uMod 


SDSS Ti-band model magnitude 




gMod 


SDSS g-band model magnitude 




rMod 


SDSS r-band model magnitude 




iMod 


SDSS i-band model magnitude 




zMod 


SDSS 2-band model magnitude 




uErr 


SDSS u-band model magnitude error 




gErr 


SDSS p-band model magnitude error 




rErr 


SDSS r-band model magnitude error 




iErr 


SDSS i-band model magnitude error 




zErr 


SDSS 2-band model magnitude error 




U R 


SDSS u — g color (corrected for ISM extinction, u — g — {uMod — 1.87rExt) — {gMod 


- 1.3SrExt)) 


gr 


SDSS g — r color (corrected for ISM extinction, g — r — {gMod — 1.38rExt) — {rMod - 


- l.OOrExt)) 


ri 


SDSS r — i color (corrected for ISM extinction, r — i — {rMod — l.OOrExt) — {iMod — 


0.76rExt)) 


iz 


SDSS i - z color (corrected for ISM extinction, iz = {iMod - 0.76rExt) - {zMod - 0.54rExt)) 


gi 


SDSS g — i color (corrected for ISM extinction, g — i — {gMod — 1.3SrExt) — {iMod — 


0.76rExt)) 


pmL c 


Proper motion in the Galactic latitude direction (in units of mas/yr) 




pmB c 


Proper motion in the Galactic longitude direction (in units of mas/yr) 




pmErr c 


Proper motion error (in units of mas/yr) 




isolated 


1 if SDSS CHILD flag is set to (isolated — — may indicate a close neighbor) 




sdss_satur 


1 if SDSS SATURATED flag is set to 1 




sdss_flags 


SDSS flags encoded as a 64-bit integer number 




theta_2mass c; 


Distance in arcscc from the 2MASS source 




J 


2MASS J-band magnitude 




JErr 


2MASS J-band magnitude error 




H 


2MASS ii-band magnitude 




HErr 


2MASS .ff-band magnitude error 




K 


2MASS K-band magnitude 




KErr 


2MASS .fC-band magnitude error 




jh 


2MASS J - H color (corrected for extinction, j — h — {J — O.Z27rExt) — {H — 0.209rExt)) 


hk 


2MASS H - K color (corrected for extinction, h — k = (H — 0.209rExt) - (K - 0.133rExt)) 


jk 


2MASS J - K color {corrected for extinction, (.7 - 0.327rExt) - (K - 0. 133rExt) ) 




a Calculatcd aftc 


r eliminating 10% of highest and 10% of lowest magnitudes. 




b This and other SDSS parameters set to NULL if no match within 2" and not in SDSS footprint, set to 


100 if in SDSS 


footprint but no match within 2" . 




c Obtained from 


thelMunn et al.l J2004h SDSS-USNOB oropcr motion catalog. 




d This and other 2MASS parameters set to NULL if no match within 2". Following: ICovev et al.1 J2007h. 


wc only match 


to 2MASS sources 


with the following 2MASS flags: rdflag — 222, blflag — 111, and ccflag — 0. 
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wavelength [nm] 

Fig. 1.— The LINEAR (dashed line) and SDSS ugriz spectral response curves (solid 
lines). The response curve for LINEAR includes the spectral quantum efficiency of the 
back-illuminated CCD, reflectivity of the aluminum mirrors, and atmospheric transmissivity 
at airmass of 1.3. The response curves for SDSS include the same atmospheric absorption. 
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Fig. 2. — The cadence of LINEAR observations measured as the time between revisits of 
the same patch of sky. From left to right, the dashed lines show 1-minute, 1-hour, 1-day, 
1-week, 1-month, and 1-year cadences. The LINEAR temporal coverage is fairly uniform, 
with a peak at ~ 1000 seconds corresponding to the main 15-min cadence. 
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Dec linear - Dec usnob [arcsec] 

Fig. 3. — A comparison of USNO-B and recalibrated LINEAR R.A. (top) and Dec single- 
epoch coordinates (bottom) for the small (solid line), LI (dotted line), and L2 setup (dashed 
line). The width of the distributions (determined from the interquartile range) puts the 
astrometric precision in LINEAR single-epoch coordinates at ~ 0.6". The systematic offsets 
in LINEAR coordinates are < 0.3" and have been subsequently corrected for. The dotted 
line at 0" was added to guide the eye. 
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Fig. 4. — The color term between the zero-point corrected LINEAR magnitude, mi nstr + Zq, 
and SDSS r-band magnitude for the small (dots), LI (solid squares), and L2 setup (open 
triangles). The solid line shows the best fit for all three setups. The rms scatter between 
the best fit for all three setups and fits to individual setups is smaller than 0.01 mag. The 
SDSS g — i color is not corrected for ISM extinction. 
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Fig. 5. — The distribution of mu near —m s d S s residuals binned in the x — y plane for the small 
(top), LI (middle), and L2 setup (bottom). The panels on the left show median values in 
each bin, and the panels on the right show rms scatter in a bin. The values are color-coded 
according to legends, with values outside the range saturating. The patterns seen on the 
left are magnitude zero-point variations that were not taken out when LINEAR images were 
flat-fielded. Therefore, the maps on the left represent "super flat-fields" that need to be 
subtracted during recalibration. 
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Fig. 6. — Similar to Fig. [5j but obtained after subtraction of "super flat-fields". The median 
values are color-coded on a smaller range than in Fig. [5] (0.02 mag vs. 0.05 mag). On average, 
the medians are now smaller than 0.005 mag, and the rms scatter is smaller than 0.05 mag. 
The higher rms scatter indicates cells with sometimes problematic readouts. 
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Fig. 7. — The medians of mn near — m S( i ss residuals in m s ^ ss bins for the small (dots), LI 
(crosses), and L2 setup (open triangles). Apart from slight non-linearity (< 0.02 mag) in the 
LI setup at bright magnitudes, recalibrated LINEAR magnitudes are linear with magnitudes 
synthesized from SDSS gri photometry. 
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Fig. 8. — The distribution of zero-point offsets between SDSS and LINEAR systems for the 
small (solid line), LI (dotted line), and L2 setup (dashed line), where the offset is measured 
for each field as the median of m Unear — m sdss residuals. The distribution is about 0.01 mag 
wide, implying that on average the LINEAR recalibration is about 0.01 mag offset from the 
SDSS photometric system. 
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Fig. 9. — The rms scatter of (mu near — m sdss ) / au near values in m sdss bins for the small (dots), 
LI (crosses), and L2 setup (open triangles). The rms scatter values greater than 1 indicate 
underestimated LINEAR magnitude errors. The pipeline-supplied errors for the small, LI, 
and L2 setup are underestimated by about 40%, 85%, and 65%, respectively. 
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Fig. 10.— The median number of observations per object as a function of equatorial 
J2000.0 right ascension and declination coordinates. The values are color-coded according 
to the legend, with values outside the range saturating. The dashed lines show ±10° of the 
ecliptic plane (solid line) and the Galactic plane is shown as a dotted line. The average 
number of observations per object within ±10° of the ecliptic plane is ~ 460, and ~ 200 
elsewhere. 
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Fig. 11. — The median number of "good" observations (observations with flags=0) as a 
function of m s d ss . The median number of "good" observations is ~ 200 for m s d ss < 18, after 
which it decreases precipitously (the 5a detection limit is at m s d ss ~ 18). About 10% of 
detections have flags> 0. 
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Fig. 12. — The cumulative distribution of x^df va l ues f° r a bright (14.5 < (mu near ) < 17, 
solid line) and faint (17 < (mu near ) < 18, dashed line) sample of LINEAR objects. The 
median of distributions is ~ 1, indicating that, on average, the uncertainties in LINEAR 
magnitudes are well determined. The high-end tail of distributions (Xpdf > 3) is populated 
by truly and spuriously variable objects (most likely dominated by the latter). 
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Fig. 13. — The median photometric error, £, as a function of magnitude. The median 
photometric error is estimated as the median of a values in a magnitude bin (dots), where 
a is the rms scatter of a time-series. The analytic description of Ti(mu near ) is given by Eq. [7] 
(solid line), and was obtained by fitting a fourth-degree polynomial to medians. The median 
photometric error increases rapidly from about 0.03 mag at 15 mag, where systematic errors 
dominate, to 0.2 mag at 18 mag, where the errors are dominated by photon statistics. 
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Fig. 14. — The completeness of the LINEAR catalog as a function of m s d ss magnitude 
(LINEAR magnitude synthesized from SDSS photometry, see Equation[6]). The completeness 
is estimated as a fraction of SDSS objects matched to LINEAR objects within 2". We find 
the object catalog to be more than 90% complete for m s dss < 19. 
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Fig. 15. — Examples of period-folded light curves from the LINEAR dataset. The left panels 
show bright ((mu near ) ~ 15) and right panels show faint objects ((mi inear ) ~ 17). The two 
top light curves are RR Lyrae stars, and other are eclipsing binary systems. The number 
in the lower left corner in each panel indicates the rms scatter around the cubic spline fit 
to the binned mean values of the light curve using 20 points (dashed line). The light curve 
features are well defined, thanks to dense LINEAR sampling, even for faint stars. 




Fig. 16. — The distribution of log(period) for ~ 7000 LINEAR variables in the SDSS u — g 
vs. g — r color-color diagram. The sample is limited to 14.5 < r < 17 and variability was 
confirmed by visual classification of light curves. Sources are binned in 0.05 mag wide color 
bins, and the median values are color-coded according to the legend on the right (values 
outside the range saturate in blue or red). Contours outline the overall distribution of SDSS 
stars (with 14.5 < r < 17) at 1%, 10%, and 50% of the peak density. 
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The LINEAR Survey Photometric Database 




Home SQL Query Object Search Catalog Matching Light Curves Documentation Contact WSG I Environment 



Now Viewing: user_sql_query 



SQL query 




Output | Online Table 



Reset 



Run Query 



Your query: "select * from source where "objectID" = 1030; " 



Matching records 

item 
1 

2 
3 
4 
5 

6 



sourcelD ra decl MJD 

1634337 115.843476 65.736482 52636.3064595 

3209021 115.843379 65.736453 52636.3498227 

6385157 115.843423 65.736468 52636.3209077 

6508400 115.843364 65.736455 52636.3353559 
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Fig. 17. — A screenshot showing the search interface for the LINEAR database. 



